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Abstract. The exponential family of random graphs is among the most widely- 
studied network models. We show that any exponential random graph model may 
alternatively be viewed as a lattice gas model with a finite Banach space norm. The 
system may then be treated by chister expansion methods from statistical nicehanics. 
In particular, we derive a convergent power series expansion for the limiting free energy 
in the case of small parameters. Since the free energy is the generating function for the 
expectations of other random variables, this characterizes the structure and behavior 
of the limiting network in this parameter region. 
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1. Introduction 

Random graphs have been widely studied (see [1, 2] for surveys of recent work) since 
the pioneering work on the independent case. The first serious attempt was made by 
Solomonoff and Rapoport [3] in the early 1950s, who proposed the "random net" model 
in their investigation into mathematical biology. A decade later, Erdos and Renyi [1] 
independently rediscovered this model and studied it exhaustively, hence the namesake 
"Erdos- Renyi random graph" . Their construction was straightforward: Take n identical 
vertices, and connect each pair by undirected edges independently with probability p. 
Many properties of this simple random graph are exactly solvable in the large n limit. 
Perhaps the most important feature is that it possesses a phase transition: From a low- 
density, low-p state in which there are few edges to a high-density, high-p state in which 
an extensive fraction of all vertices are joined together in a single giant component. 

The Erdos-Renyi random graph, while illuminating, is a poor model for most real- 
world networks, as has been argued by many authors El El, and so it has been 
extended in a variety of ways. To address its unreahstic degree distribution, generalized 
random graph models such as the configuration model [H] and the multipartite graph 
model [9] have been developed. However, these models have a serious shortcoming, in 
that they fail to capture the common phenomenon of transitivity exhibited in social and 
biological networks of various kinds. 

The main hope for progress in this direction seems to lie in formulating a model that 
incorporates graph structure in more detail. A top candidate is the exponential random 
graph model, in which dependence between the random edges is defined through certain 
finite subgraphs, in imitation of the use of potential energy to provide dependence 
between particle states in a grand canonical ensemble of statistical physics. These 
exponential models were first studied by Holland and Leinhardt [10] in the directed case, 
and later developed extensively by Frank and Strauss [H], who related the random graph 
edges to a Markov random field. More developments are summarized in [121 UHl E] • 

The past few years have witnessed a surge of interest in the study of the limiting 
behavior of exponential random graphs. A major problem in this field is the evaluation 
of the free energy, a quantity that is crucial for carrying out maximum likelihood and 
Bayesian inference. A particular motivation for people in the statistical mechanics 
community to study the free energy is that it provides information on phase transition 
in these sophisticated models. Many people have made substantial contributions in 
this area: Haggstrom and Jonasson [15] examined the phase structure and percolation 
phenomenon of the random triangle model. Park and Newman fWi [T7] constructed 
mean-field approximations and analyzed the phase diagram for the edge-two-star and 
edge-triangle models. Borgs et al. [18] established a lower bound on the largest 
component above the critical threshold for random subgraphs of the n-cube. Bollobas et 
al. [19] showed that for inhomogeneous random graphs with (conditional) independence 
between the edges, the critical point of the phase transition and the size of the giant 
component above the transition could be determined under one very weak assumption. 
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Bhamidi et al. [20j focused on the mixing time of the Glauber dynamics and proposed 
that in the high temperature regime the exponential random graph is not appreciably 
different from the Erdos-Renyi random graph. Dembo and Montanari [21] discovered 
that for the Ising model on a sparse graph phase transitions and coexistence phenomena 
are related to Gibbs measures on infinite trees. Using the emerging tools of graph limits 
as developed by Lovasz and coworkers [22], Chatterjee and Diaconis [23] gave the first 
rigorous proof of singular behavior in the edge-triangle model. They also suggested 
that, quite generally, models with repulsion exhibit a transition qualitatively like the 
solid/fluid transition, in which one phase has nontrivial structure, as distinguished from 
the "disordered" Erdos-Renyi graphs. Radin and Yin [21] derived the full phase diagram 
for a large family of 2-parameter exponential random graph models with attraction and 
showed that they all contain a first order transition curve ending in a second order 
critical point (qualitatively similar to the gas/liquid transition in equilibrium materials). 
Aristoff and Radin [25] considered random graph models with repulsion and proved that 
the region of parameter space corresponding to multipartite structure is separated by a 
phase transition from the region of disordered graphs (proof recently improved by Yin 

As is usual in statistical mechanics, we work with a finite probability space, and 
interpret our results in some more sophisticated limiting sense. We consider general 
fc-parameter families of exponential random graphs. 

• A complete graph Kn on n vertices consists of a vertex set (|V^| = n) and an 
edge set En (|-E„| = (2))- A vertex pair e is a two-element subset of Vn, and the 
set of all vertex pairs constitute En- 

• Qn is the set of simple graphs G on n vertices, where a graph G with vertex set 
V{G) = is simple if its edge set E{G) is a subset of En- 

• ifi, ...,Hk are pre-chosen finite simple graphs. Each Hi has rrii vertices (2 < rrii < 
m) and Pi edges {1 < Pi < p) ■ In particular. Hi is K2 (i.e. a single edge). 

• A vertex map / : V{Hi) y{G) is a homomorphism if the induced edge map 

: E{Hi) — En sends E{Hi) into E{G). t{Hi,G) is the density of graph 
homomorphisms H^ — )■ G: 

_ |hom(g„G')| 

t{Hi,G) - — , (Ij 

\hom{Hi, Kn)\ 

where the denominator \hom{Hi, Kn)\ = n™' counts the total number of mappings 
from V{Hi) to V{G). 

• (3 = (/3i, ...,/3fc) are k real parameters. They tune the influence of the pre-chosen 
graphs Hi,...,Hk. 

Let T^{G) be the weighted sum of graph homomorphism densities t{Hi, G): 

k 

TP{G) = Y,P,t{H,,G). (2) 

i=l 
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There is a probability mass function that assigns to every G 

F^^iG)=exp{n\T^iG)-rn)), (3) 
where tp^ is the normahzation constant, i.e., it satisfies 

5^exp(n2T^(G))=exp(nV9- (4) 

The rest of this paper is organized as follows. In Section [2] we show that the 
general exponential random graph model may alternatively be viewed as a lattice gas 



model with a finite Banach space norm (Propositions 2.3 and 2.4). This transforms the 
probability model into a statistical mechanics model (Theorem 2.5). In Section [3] we 
apply cluster expansion techniques [27] and derive a convergent power series expansion 
(high-temperature expansion) for the limiting free energy in the case of small parameters 



(Theorem 3.5). Finally, Section |4] is devoted to concluding remarks. 



(Jp. 



2. Alternative view 

In this section we will transform the exponential random graph model into a lattice gas 
model with a finite Banach space norm. We begin by presenting an alternative view of 
the homomorphism density t(ifi, G) ([!]). 

Definition 2.1. Let G G Let a be the indicator function of E{G). For every vertex 
pair e = {i,j}, 

1, an edge exists between vertices i and j; 
0, otherwise. 

Definition 2.2. Let X C En- Fix a finite simple graph H with p edges. Define the 
exact graph homomorphism density d{H, X) by 

\hom{H,Kn)\ 

where the numerator \ehom{H, X)\ counts the number of homomorphisms / : V{H) — )■ 
Vn whose induced map : E{H) — )• En sends E{H) onto X. It is clear that d{H,X) is 
finite-body: d{H,X) = for |X| > p. 

Proposition 2.3. Let G G Qn- Let a be the indicator function of E{G). Fix a 
finite simple graph H. The graph homomorphism density t{H, G) has a lattice gas 
representation 

t{H,G) = J2 d{H,X)ax, (6) 

X<ZE„ 

where ax = Ileex ^e- 

Proof. By Definition |2.2| the graph homomorphism density ([T]) is given by 

t{H,G)= J2 d{H,X). (7) 

XCE(G) 

Our claim easily follows once we realize that = 1 if and only if X C E{G). □ 




Figure 1. Constructing a lattice gas representation for the graph homomorphism 
density t{H, G). 



The above construction might be better explained with a concrete example. Take 
H a two-star (with m = 3 vertices and p = 2 edges). Let G (see Figure [T]) be a finite 
simple graph on 4 vertices (n = 4): A, B, C, and D. There are 18 homomorphisms and 
43 = 64 total mappings from V{H) to V{G). By ([T]), the homomorphism density of H 
in G is t{H,G) = 18/64. 

The density t{H, G) may be derived through a lattice gas representation as well. 
For notational convenience, we denote d{H,X) by J{X). The image of if in G under a 
homomorphic mapping is either an edge or a two-star. The exact homomorphism density 
for an edge is 2/64, and we have J{{A,B}) = J{{A,G}) = J{{A,D}) = J{{B,G}) = 
J{{B,D}) = J{{C,D}) = 2/64. The exact homomorphism density for a two-star is 
also 2/64, and we have J{{A, B}, {B, C}) = J{{B, A}, {A, G}) = J{{A, G}, {G, B}) = 
J{{A,B},{B,D}) = J{{B,A},{A,D}) = J{{A, D}, {D, B}) = J{{A,G}, {G, D}) = 
J{{G,A},{A,D}) = J{{A,D},{D,G}) = J{{B,G}, {G, D}) = J{{G,B},{B,D}) = 
J{{B,D},{D,G}) = 2/64. For all other X C E^, we have J(X) = 0. The indicator 
function cr of this particular E{G) is given by: cr{A,B} = '^{a,d} = (^{3,0} = '^{s.-D} = 1 
and cr^A,c} = c"{c,d} = 0. Therefore the valid images of if in G are: edges AB, AD, BG, 
BD, each carrying density 2/64; two-stars {AB, BG), {AB, BD), {GB, BD), {BA, AD), 
{AD, DB), each carrying density 2/64, making the combined density YIixclea J{^)^x = 
18/64. 

As n gets large, it would seem hard to keep track of all the densities d{H, X) in 
the lattice gas representation, nevertheless, we will show that the pinned densities have 
a universal upper bound. 

Proposition 2.4. Let G E Qn- Fix a finite simple graph H with m vertices. Fix a 
vertex pair e = of G. Denote by te{H,G) the sum of the exact homomorphism 

densities d{H,X) with e G X C E{G). Then in the lattice gas representation (cf. 
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Proposition 2.3), we have 



t,iH,G)= Yl diH,X)<'^^^^^. (8) 

X:e&XCE{G) 

Proof. The homomorphisms under consideration must satisfy that the image of V{H) 
in V{G) contain vertices i and j of G. To count these homomorphisms, we regard 
such a mapping as consisting of two steps. Step 1: We construct vertex maps from 
V{H) to V{G): First select two vertices of H to map onto i and j, of which there are 
m(m — 1) choices. Then map the remaining vertices of H onto G, of which there are 
n™~^ ways. Step 2: We check whether these vertex maps are vahd homomorphisms (i.e. 
edge-preserving). The number of homomorphisms is thus bounded by m(m — l)n'^~'^. 
Our claim easily follows once we recall that the total number of mappings from V{H) 
to V{G) is n'". □ 

Our next theorem formulates the exponential random graph model as a lattice gas 
model. The lattice is En, the set of all vertex pairs e of Vn- For each site e = {i,j}, 
we attach a lattice gas variable cTg which takes on the value 1 or 0. This specifies a 



simple graph G on Vn by Definition 2.1 The Hamiltonian H{a) is a weighted sum of 



graph homomorphism densities t(Hi, G) and varies according to the structure of the 



graph G (cf. Proposition 2.3). It has a finite Banach space norm which depends on 



the universal upper bound for the pinned densities d{Hi,X) (cf. Proposition 2.4). This 
model is somewhat unconventional in the sense that the underlying lattice is "infinite- 
dimensional" (a vertex pair ci is a nearest neighbor of another vertex pair 62 as long as 
they have a vertex in common, so the number of nearest neighbors of a given site grows 
with n, rather than being fixed at 2'^ in a rf-dimensional lattice), yet the associated 
interaction is finite-body. A key interest in this model is its behavior in the large n 
limit. For that purpose, we will pay special attention to the limiting free energy as it is 
the generating function for the expectations of other random variables. 

Theorem 2.5. The general exponential random graph model ^ is equivalent to a 
lattice gas model with a finite Banach space norm. 

Proof. The Hamiltonian is the negative of the exponent of the probability mass function 
(|3| (without normalization): 

k 

/f(a) = -n2^At(if„G'), (9) 
1=1 



where a is the indicator function of E{G). By Proposition 2.3, it has a lattice gas 
representation 

k 

H{a) = -n^J2f^^ Yl diHi,X)ax. (10) 

i=l XCEn 
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k 

K{X)=n''Y,P'AH,,X). (11) 



i=l 



The Hamiltonian notation may then be simphfied, 

H{a) = -J2 K{X)ax. (12) 

The interactions K = {K{X) : X C En} vary with the tuning parameters Pi and form 
a Banach space B. Define the i3-norm of K by 

||ir||=sup J2 I^WI- (13) 

^ X-.e&X 



By Proposition 2.4 



\K\\ < m{m - (14) 



This Banach space construction will be useful later on. □ 

Notice that K{X), being a weighted sum of d{Hi,X), satisfies the finite-body 
property: K{X) = for |X| > p (its importance to be illustrated in the proof of Lemma 
3.6). Moreover, summing over G G (|4]) is equivalent to summing over indicator 



functions a of E{G) C En- 

J2 exp (n^ Yl ^) ) = E (-^(^)) • (15) 

Define ip'^ to be the limiting free energy of the random graph model, i.e., ip'^ = 



lim^^oo V-f. By (|4|, (p|, and (pl, we have 

= hm ^log I > exp I > ] K{X)ax 1 1 • (16) 



We normalize the sum over a for the ease of cluster expansion. Henceforth 

will denote the normalized sum and satisfy 1 = 1, i.e., = 2~l-^"l = 

2~^^ 'Ylia- Define the partition function W by 

norm norm / \ 

= ^ exp {-H{a)) = ^ exp ^(^)^x ■ (17) 

a- a \xCEr, J 

According to standard statistical mechanics, the limiting free energy 0^ of the lattice 
gas model is then given by 

/ = hm -^logW. (18) 

n— s>oo \En\ 
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We explore the relationship between the two limiting free energies ip^ and (j)^: 

= lim ^ log (2('^)w] = ^ I log2 + lim 7^ \ogw] (19) 

n-5-oo V / 2 \ n-5-oo (^j / 

= I flog 2+ lim ^logl^) = J(log2 + 0^). (20) 

The two interpretations of the limiting free energy are thus not appreciably different, 
and we may interpret it in either way to help with the understanding of the structure 
and behavior of the limiting network. 

3. Cluster expansion 

In this section we will apply cluster expansion techniques to derive a convergent power 
series expansion (high-temperature expansion) for the limiting free energy in the case 
of small parameters. The cluster expansion expressions presented here are completely 
rigorous for finite models, and may be interpreted in some more sophisticated limiting 
sense. We begin by introducing some combinatorial concepts. A hypergraph is a set of 
sites together with a collection F of nonempty subsets. Such a nonempty set is referred 
to as a hyper-edge or link. Two links are connected if they overlap. The support of 
a hypergraph is the set UF of sites that belong to some set in F. A hypergraph F is 
connected if the support of F is nonempty and cannot be partitioned into nonempty 
sets with no connected links. We use Fc to indicate connectivity of the hypergraph 
Fc. Our first proposition gives a cluster representation for the partition function of the 
exponential random graph model. 

Proposition 3.1. Let W be the partition function of the exponential random graph 
model on n vertices ( [l7| ). Then W has a formal cluster representation 

^ = E n (21) 

A AfGA 

where: 

• A is a set of disjoint subsets N of En- 

• = Eur.=^EH7nxer.(e^^^^'^^ - !)• 

Proof. We rewrite exp (XlxcSn -^(^)'^x) as a perturbation around zero interaction, 

norm norm 

^ = E n (1 + - 1) = E E n (^''^''^''" - 1) ' (22) 

a XCE„ a V X& 

where F is a set of subsets X of i?„. 



We are going to organize the sum over hypergraphs in (22) in the following way. 
Let be a possible support for a connected hypergraph. Let A be a disjoint set of such 
sets A^. Let S" be a function that takes A^ G A to a hypergraph with support A^, i.e., 
US'(A^) = A^. Then summing over hypergraphs F is equivalent to summing over A and 
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functions 5* with the appropriate property. Furthermore, the product over in A and 
the hnks in ^(A^) is equivalent to the product over the corresponding F. We have 

norm 

'•'^EEEII n (e™"-l). (23) 

<T A S NeAXeS(N) 

By independence, the sum over a can be factored over A. Denote by Yl'a°\iT 
normahzed sum restricted to subset A^, i.e., ~ 2^1^' J2cr\N- Similarly, denote 

by Yla°\e^ normalized sum restriced to vertex pair e. This gives 

(norm \ norm 
IlEl EIIE n (e"™--!). (24) 
e^UA a\e I S NeA a\N X€S{N) 



Because of the normalization, Xlafe™ 1 = (24) can be simplified 



norm 



'^'^EEIIE n (25) 

A 5 NeA (7\N X(^S(N) 

Rearranging the terms by the distributive law, we have 

norm 

= E n E E n i-"""'" - 1) ■ (26) 

A NeAur^=N a\N XgFc 

Our claim follows once we recall the definition of w^. □ 
Notice that ([2T|) has a graphical representation: 



E n(i-^(^-^^))^^i---^^" (27) 

n=0 ■ Ni,...,Nn {i,j} 
oo ^ 

E- E c{N,,...,N^)wn,---wn„, (28) 



n=0 Ni,...,N„ 



where 



c(iVi,...,iV„) = ^ II (-t(iV„iV,)), (29) 

R {i,j}&R 

i? is a graph with vertex set {1, n}, and 

t{N.,N,) = i ^ if and AT, overlap; 
10 otherwise. 

This alternate expression of the partition function W facilitates the application of 
cluster expansion ideas, which roughly summarized, state that a sum over arbitrary 
graphs can be written as the exponential of a sum over connected graphs. Taking the 
logarithm of the partition function (28) thus replaces the sum over graphs by the sum 
over connected graphs. The log operation is physically significant in that the resulting 
connected function logVF is proportional to the limiting free energy 0^ (18). A detailed 
explanation of this phenomenon may be found, for instance, in a survey article by Faris 
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Proposition 3.2. Let W be the partition function of the exponential random graph 
model on n vertices (28). Then the connected function logiy is given by the cluster 
expansion 

logW = J2- E C{N,,...,N^)wr,,---wr,„, (31) 



w. 

n=l Ni,...,N, 



where 



CiN,,...,N^) = J2 n i-m,Nj)), (32) 
Rc {i,j}&Rc 

and Rc is a connected graph with vertex set {1, ...,n}. 

Now that we have derived an explicit expression for the connected function log W, 
we explore criteria that guarantee the convergence of this expansion. This provides 
information on the limiting free energy (p^ (18) and characterizes the structure and 
behavior of the limiting network. The celebrated theorem of Kotecky and Preiss says 
that if the interaction is sufficiently weak, then the cluster expansion for the pinned 
connected function converges. 

Theorem 3.3 (Kotecky-Preiss [29]). Consider arbitrary family of activities > 0. 
Suppose there are finite -Bat > such that for all A^'o, 

J2tiN,No)vNBN<\ogBNo. (33) 

N 

Then the pinned connected function has a convergent power series expansion in the 
polydisc \wn\ < vn- 



oo ^ 

Yl \C{Ni,...,Nn)\\wNA---\wNj<VNBM. (34) 

n=l ' Ni,...,N„:3iN,=N 

Remark. In application it is convenient to take = M'^' with M > 1. With this 
choice of B^ the Kotecky-Preiss condition is equivalent to the condition that (33) holds 
for all one point sets Nq = {e}: 

J2 wtvMI^I < logM. (35) 

N-.e&N 

This reduced version of the Kotecky-Preiss condition will be used throughout the rest 
of this section. 



At first sight, the Kotecky-Preiss condition (35) seems very abstract and difficult 
to verify. It is a weighted sum of activities v^- pinned at the vertex pair e, and each 
is an upper bound for \wn\, whose expression involves a hypergraph decomposition and 
is rather complicated by itself (cf. Proposition 3.1). The following proposition gives a 
handy criterion for weak interaction in the small parameter region (^^^^ |A| small). 



Proposition 3.4. Fix M > 1. Take wjy as in Proposition 3.1 It is clear that 

urc=N X&Fc 
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Consider the interaction K with the Banach space norm \\K\\ (cf. Theorem 2.5). 
Suppose Yl'i=i lAI is smaU: 

\\K\\ < ,n(,n - i)±m < .,.J°^':''r',r.-,.. - (37) 



1=1 



2(Mp)p(l + {p- l)logM)' 



Then (35) holds for every vertex pair e. 



Remark. The maximal region of parameters {Pi} is obtained by setting 
logM : 



-p + 5p2 — Ap 



(38) 



2p{p - 1) 

Our ultimate goal is to examine convergence of the limiting free energy 0^ (18), 
which is proportional to the connected function logly (31). As the Kotecky-Preiss 
result concerns convergence of the pinned connected function (34), it appears to be 
inapplicable. However, our next theorem shows that pinning is actually the central 
ingredient that ties these two seemingly unrelated issues together. 

Theorem 3.5 (Main Theorem). Fix M > 1. Consider the interaction K with the 
Banach space norm (cf. Theorem 2.5). Suppose Yli=i lAI small (37). Then the 
limiting free energy 0^ ( 18 ) is analytic in P and the rate of convergence is uniform. 



Proof. By Proposition 3.4, (35) holds for every vertex pair e. The pinned connected 
function (34) thus converges absolutely by Theorem 3.3, which further implies 

oo ^ 



|C(iVi,...,iV„)|KJ---kivJ (39) 

N<ZE„n=l '"' Ni,...,N„:3iN,=N 

< J2 ^ivMl^l <J2J2 ^ivMl^l < iKllogM. (40) 

We conclude that the limiting free energy 0^ is absolutely convergent and bounded 
above by log M. □ 



The rest of this section is devoted to the proof of Proposition |3.4[ The weighted 
activity sum in the Kotecky-Preiss weak interaction condition (35) is rewritten as a 
power series, whose terms are then shown to be exponentially small under (37) by a 
series of lemmas. 



Proof of Proposition 3.4 We notice that when \\K\\ is small (say \\K\\ < |), 
q\k(x)\ _ ^ < 2\K{X)\ by the mean value theorem. For UP^ = A^, \N\ < ^ |X| with X 
in Pc. We have 



N:eeN 



< 



(41) 
(42) 



TeieeurcXerc 
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We say that a hypergraph Tc is rooted at the vertex pair e if e G UFc. Let a„(e) be the 
contribution of all connected hypergraphs with n links that are rooted at e, 

an{e)= n2|^WI^""- (43) 

eeurc:|rc|=nXerc 



Then 



N:eeN n=l "^^-^^ n=l 

It seems that once we show that a„ is exponentially small, the power series above will 
converge, and our claim might follow. To estimate a„, we relate to some standard 
combinatorial facts [30] . □ 

Lemma 3.6. Let a„ be the supremum over e of the contribution of connected 
hypergraphs with n links that are rooted at e. Then a„ satisfies the recursive bound 

an<2\\K\\MPY.(t) E «n.---a., (45) 

A;=0 ^ ^ a„j,...,a„^:ni+...+ns,+l=n 

for n > 1, where (^) is the binomial coefficient. 

Proof. We linearly order the vertex pairs e in En and also linearly order the subsets X of 
En- For a fixed but arbitrarily chosen e in we examine (43). Write Fc = {^i} U Fj, 
where Xi is the least X in Fc with e G Xi. As -ft'(Xi) 7^ only for \Xi\ < p, 

a„(e) < 2\\K\\MPJ2 H 2|i^(X)|Ml^l. (46) 

The remaining hypergraph Fj has n — 1 subsets and breaks into k : k < p connected 
components (which again follows from the finite-body property of K). Say they are 
Fi, Ffc of sizes ni, Uk, with ni + ... + rik = n — 1. For each component Fj, there is a 
least vertex pair through which it is connected to Xi, i.e., Cj G Xi is a root of Fj. As 
1-^1 1 < P and Fj consists of k components, the number of possible choices for the root 
locations is at most (^) . We have 

anie)<2\\K\\Mpf2(t) Yl a„,---a„,. (47) 

k=0 ^ ^ ani,..-,ani.-ni+...+nk+l=n 

Our inductive claim follows by taking the supremum over all e in En- Finally, we look 
at the base step: n = 1. In this simple reasoned above, we have 

ai = sup 

eeurc:|rc|=ixerc 

and this verifies our claim. □ 
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Clearly, J^N-eeN vnM^^^ will be bounded above by J2'^=i '^n, if 



I — n V / 



fc=0 ^ ^ a„j,...,a„j,:ni+...+rafe+l=n 

for n > 1, i.e., equality is obtained in the above lemma. Observe that when n = 1, the 
empty product (corresponding to = 0) is the only nonzero term on the right side of 
(49), so we have di = 2\\K\\Mp, matching the bound in (48). 

Lemma 3.7. Consider the coefficients a„ that bound the contributions of connected 
and rooted hypergraphs with n links. Let w = Yl'^=i (^nZ"" be the generating function 
of these coefficients. The recursion relation ( [49] ) for the coefficients is equivalent to the 
formal power series generating function identity 

w = 2\\K\\MPz{l + wf. (50) 



Proof. Notice that {1 + w)'p = Y1=q thus 

V 

w = 2\\K\\MPzJ2 



k=0 



Writing completely in terms of z, we have 

p / 

P 
k 



Y,anz'' = 2\\K\\MPY,[ 

n=l k=0 ^ 



.,a„. :rai+...+nfc+l=n 



(51) 



(52) 



We compare the coefficient of z on both sides of (52): On the left it is given by di, and 
on the right it is given by 2||iir||M^ times the empty product, thus di = 2\\K\\Mp. Our 
general claim follows from term-by-term comparison. □ 



Lemma 3.8. If w is given as a function of 2; as a formal power series by the generating 
functi 

\z\ < 



function identity (50), then this power series has a nonzero radius of convergence 



2||X||(Mp)p- 



Proof. Without loss of generality, assume ^ > 0. Set zi = 2\\K\\MPz. Solving (50) for 
Zi gives Zi = w/{l + wy. As Zi goes from to (p — ly^^ /pP, the w values range from 
Otol/(p-l). □ 



Proof of Proposition \3.4\ continued. We notice that in the above lemma, w = 
Yl'^=i = l/(p ^ 1) corresponds to 2117^1171^*^2; = {p — iy~^/pP, which implies that 
for each n, 

d^ < {2\\K\\{MpYr {p - ly^'^^^-'^'^K 

Gathering all the information we have obtained so far. 



t;;vMl^l< 5^(2||ir||(Mp)T(p-l)" 



-(l+(p-l)n) 



N:eeN 



n=l 



2\\K\\iMp)P 

(p-l)P 
_ 2||A-||(A/p) 



- <logM. 



by dSTl). 



(53) 
(54) 

(55) 

□ 
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4. Concluding remeirks 

This paper reveals a deep connection between random graphs and lattice gas (Ising) 
systems, making the exponential random graph model treatable by cluster expansion 
techniques from statistical mechanics. We show that any exponential random graph 
model may alternatively be viewed as a lattice gas model with a finite Banach space 
norm and derive a convergent power series expansion (high-temperature expansion) for 
the limiting free energy in the case of small parameters. Since the free energy is the 
generating function for the expectations of other random variables, this characterizes the 
structure and behavior of the limiting network in this parameter region. We hope this 
rigorous expansion will provide insight into the limiting structure of exponential random 
graphs in other parameter regions and shed light on the application of renormalization 
group ideas to these models. 
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